Dynamical Janssen effect on granular packings with moving walls 
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Apparent mass (M app ) measurements at the bottom of granular packings inside a vertical tube in 
relative motion are reported. They demonstrate that Janssen's model is valid over a broad range of 
velocities v. The variability of the measurements is lower than for static packings and the theoretical 
exponential increase of M app with the height of the packing is precisely followed (the corresponding 
characteristic screening length is of the order of the tube diameter). The limiting apparent mass at 
large heights is independent of v and significantly lower than the static value. 

PACS numbers: 45.70.-n, 81.05.Rm, 83.80.Fg 
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Dense, dry vertical particle flows in channels are fre- 
quently encountered in industrial processes such as the 
emptying of silos or pneumatic transport y]. Particle 
fraction and velocity variations in these flows depend 
largely on interaction forces either between the parti- 
cles or between the particles and the channel walls [2, • 
The present work studies experimentally such forces in 
the case of a granular packing globally at rest but in rel- 
ative motion with respect to an upwards moving tube. 
These force distributions have been extensively stud- 
ied by many authors for static or quasi-static packing 
in tubes ff E l 0. 0. or other configurations 
[ill H2I Il3l Il4|. The pioneering work of Janssen demon- 
strated in particular that the effective weight at the bot- 
tom end of a container reaches exponentially a limit when 
the height of a static packing increases: this is due to the 
shielding effect of contact forces between grains redirect- 
ing the weight towards the side walls. Vanel et al. 01 
measured the apparent weight of packed beads moving 
very slowly down a vertical tube (v = 20^ms -1 ): these 
experiments indicate that the weight increases during the 
displacement before reaching a constant value but the 
influence of the height of the packing was not investi- 
gated. In this work we demonstrate experimentally that 
Janssen's results can be generalized to a grain packing 
moving with respect to solid walls at velocities up to sev- 
eral cms -1 . 

A packing of 2 mm diameter glass beads (density 
p = 2.53 10 3 kgm -3 ) is contained inside a vertical glass 
tube of length 400 mm and internal diameter D = 30 mm. 
This tube can be moved up and down at a chosen veloc- 
ity v ranging from 10/ims -1 to 35mms -1 . A cylindrical 
piston of height 80 mm and diameter 29 mm is inserted in 
the bottom part of the tube: it is carefully aligned with 
it and does not touch the walls during the motion. The 
bottom end of the piston is screwed onto a strain gauge 
sensor connected to a lock-in detector allowing one to 
determine the force of the beads on the upper horizon- 
tal surface of the piston (Fig. This value is divided 
by the gravitational acceleration g to obtain an apparent 
mass -Mapp determined with an uncertainty of ±0.1g. 
This apparent mass can be compared to the actual mass 
M of the grains measured by electronic scales prior to 



the experiment, with a precision of ± 0.01 g. In addition, 
qualitative observations are realized using a video camera 
fitted with a macro lens and connected to a video tape 
recorder. The experiments were performed at a relative 
humidity H = (50 ± 5)%. For each experimental run, the 




FIG. 1: Sketch of the experiment. 

tube is first brought to its low position and filled with a 
given mass of beads M, using a funnel placed at the top of 
the tube. The relative motion between the grains and the 
walls is then induced by moving the tube upwards over 
a distance Az ~ 70 mm. In the initial and final phases 
of the motion, the acceleration is equal to 0.04 ms -2 (or 
to 0.055 ms -2 for v = 35mms -1 ) and the velocity v is 
constant in between. At the fastest velocity, the acceler- 
ation and deceleration distances are of order 11mm and 
decrease to less than 10 -2 mm for v <lmms . One 
expects that, in the constant velocity phase, the relative 
motion of the grains with respect to the walls (and the 
force distribution) will be identical to those for a steady 
downwards grain flow in a tube at rest. The two con- 
trol parameters of the experiment are the mass M of 
the grains poured into the tube and the tube velocity v. 
Compared to the flow of grains inside a static tube, the 
relative motion of grains and of the walls is the same but 
both the accelerations in the initial and final phases and 
the flow of air differ. Air flow should however have a 
weak influence since pressure gradients in these station- 
ary compact grain flows are generally small. 

Figure [21 displays variations of the apparent mass 
M app as a function of time for seven identical exper- 



2 



iments. The mass of the grains initially poured into 
the tube is M = 300 g and corresponds to a height 
h = (262 ±1) mm of the grain packing. Values of the ap- 
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FIG. 2: Variation of the apparent mass M apP of a grain pack- 
ing of total mass M = 300 g as a function of time: ( 1) right af- 
ter filling the tube, (2) tube moving upwards (» = 4mms~'), 
(3) tube stopped. Various line patterns correspond to differ- 
ent experiments realized with same control parameters. 
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FIG. 3: Apparent mass as a function of total mass 
M of grain packing for different tube velocities: (0,^) 
ii = 0.02mms _1 , (□,■) u = 0.2mms _1 , (o, •) » = 2mms" 1 1 
(A, A) « = 20mms _1 , (V, T) D = 34mms _1 . Each line corre- 
sponds to an exponential fit of the data points at each ve- 
locity. Dark symbols: limit value at the end of region 
(2) (moving tube) (AMfpp 1 = ±lg). Open symbols: appar- 
ent mass M| P p' in region (3) (stopped tube) (AM'" = ± 8 g). 
Dashed thick line: M aDD = M. 



parent mass M app of the static packing right after filling 
the tube (region (1) in Fig. are dispersed: the devia- 
tion is of the order of ± 25% of the mean value for large 
heights. Similar observations are reported by other au- 
thors 0, 0] : they are explained by the heterogeneous 
stress network created when pouring the grains into the 
tube. Arches of grains appear at random inside the pack- 
ing and screen the weight of other grains located above. 
These effects vary strongly from one sample to another, 
resulting in large variations of Af app . 

At the onset of the tube motion, the apparent mass 
decreases abruptly and a random motion of the beads is 
observed. Then, M app increases slightly during the con- 
stant velocity phase (region (2) in Fig. [5J and reaches a 
limiting value referred to in the following as M^£. A 
striking feature is the fact that does not vary by 

more than ± 5% of its mean value from one sample to 
another. Figure El also shows that the time dependence 
is the same for all experiments within experimental un- 
certainty. When the tube stops (region (3) in Fig. |2), 
the apparent mass increases sharply to a constant value 
M| P p* that is still much lower than in phase (1). The rel- 
ative dispersion of these last static values is again quite 
large (± 15%). Moreover, one observes that the values of 
the apparent mass in region (1) and region (3) are not 
correlated. These experiments were repeated for different 
total masses M of grains and different tube velocities v. 
The respective experimental mean final values of M app 
after it has become constant in region (2) and after the 
tube has stopped (region 3) are plotted in Fig. El 

The apparent mass M^g" at the end of the moving 
phase (dark symbols), is almost independent of the ve- 



locity for all total masses M of grains in the tube. At 
low values of M, one has ~ M (there is no screen- 

ing effect). Then the curve levels off and reaches a 

limit value = (34.5 ± 1) g, independent of v. 

After the tube motion has stopped (region 3), the 
apparent mass increases by 20 to 60% (open symbols) 
as already noted on Fig. |21 This increase is larger at 
higher velocities, but the final value is always lower than 
in region (1). The global shape of the variation of M| Pp * 
with M is qualitatively similar to that obtained for M^pp , 
although the dispersion of the data points is much higher. 

Video recordings of the grain packing during the ex- 
periments allowed one to estimate variations of the mean 
global particle fraction c in the sample: no systematic 
motion of the upper interface was detected and c re- 
mained constant with c = (64 ± 0.5)% in all experiments. 
These results indicate that the motion of the walls re- 
organizes the internal structure of the column: it allows 
to reach always a same statistical dynamical equilibrium 
in spite of the variability of the force distribution in the 
initial grain packing. This reorganization results from 
very minute motions since no variation of c is detected. 
Bringing the walls to rest perturbs the stress distribu- 
tion and leaves the system "frozen" in a state which may 
vary significantly from one sample to another. As could 
be expected, the perturbation is larger at high velocities. 

The curves of Fig. |3| follow qualitatively the predic- 
tions of Janssen's model for static packing. We investi- 
gated therefore whether the quantitative predictions were 
also verified. The model assumes that the weight of the 
beads is balanced by variations of the stress a zz in the 
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vertical direction z and by wall friction with: 

da zz 



dz 
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peg, 
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where c is the local particle fraction of beads and 
A = D/AfiK a characteristic length. The coefficient K 
characterizes the deflection of vertical stresses into hori- 
zontal radial ones in the packing and [i is the Coulomb 
friction coefficient. Integrating between the bottom 
(z — zq) and the surface (z = 0) of the packing where 
c 2Z (0) = 0, gives the stress a zz (zo) on the piston: 



Since M app : 
lows that: 



<?zM =pcgX(l-e- z °/ x ). (2) 
: cr zz (z )iTD 2 /4g and M = pcz nD 2 /4, it fol- 
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where M m is the predicted asymptotic value of the ap- 
parent mass M ap p when M is increased: 



Moo = pc\irD 2 /A. 



(4) 



The solid and dashed lines in Fig. |31 represent respec- 
tively the best fit with Eq. @ of experimental data points 
obtained during and after the motion of the tube. The 
fits are excellent in the first case and fair in the second 
one due to the larger dispersion of the data points. These 
fits allow one to determine the asymptotic apparent mass 
Afoo an d the characteristic length A of Janssen's model. 
Variations with the tube velocity v of A values corre- 
sponding to both sets of data are displayed in Fig. In 



50 



45 



E 40 

E. 



35 



30 



25 - 



D 


= 30 mm 










d 


= 2 mm 


















^1 








1 


1 ^ " " 1 


1 




1 


1 


1 








i 


■ 


• • 





0.01 



I 

0.1 



- 1 1 — I I I I 1 1 1 



1 10 

V (mm s" 1 ) 



I 

LOO 



FIG. 4: Janssen's length variations with the tube velocity: 
(T) during the flow, (■) after the tube stops. 

the moving phase, A remains at the nearly constant value 
A = (30 ± 2) mm. This value is close to the tube diame- 
ter (D = 30 mm) and two or three times lower than esti- 
mates in region (1) from Eq. (0}. After the tube has been 
stopped (region 3), A is nearly the same as in region (2) 



at the lowest velocity but increases by 30% at the high- 
est one. This confirms that perturbations of the stress 
distribution induced by stopping the tube are larger at 
higher velocities. A part of the arches and of the con- 
tacts between grains probably disappear, which reduces 
the value of the coefficient K (and therefore increases A). 

Another important issue is the dynamics of the grain 
rearrangement while the tube is moving. This can be 
inferred from variations of the apparent weight during 
the motion of the tube. Figure displays superimposed 
variations of M app /M^ as a function of the normal- 
ized displacement z/D for several velocities ranging from 
40/ims -1 to 34mms _1 . The total mass of grains is equal 
to M = 300 g for which = = 34.5 g. A common 
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FIG 



Variation of the normalized apparent mass 
,/Mfpp with the normalized displacement z/D, at ve- 



locities 40 fim s < v < 34 mm s 

increasing trend is followed and is well adjusted by an 
exponential variation (dashed line): 



M„. 



M a d p y p n 



1 — A exp 



az 



(5) 



with A = 0.48 and a = 2±0.2. The characteristic relax- 
ation length towards the limit value in Fig. [S] is 
equal to D/2. The sharp drop of M app right after the on- 
set of the motion implies that grains are initially dragged 
upwards (a small transient motion is observed on video 
recordings); then, the packing rearranges and more con- 
tact paths build-up towards the sensor at the bottom. 

The influence of the tube diameter D has been inves- 
tigated by performing the experiments on a tube of larger 
diameter D = 40 mm (beads of diameter d = 3 mm are 
used to keep D/d roughly constant). As for D = 30 mm, 
the variation of M app /Al^ with z/D is almost indepen- 
dent on the tube velocity (Fig. 0) . M app increases also 
exponentially with the displacement z, following Eq. (JSJ 
with the same value of the parameter a, as for D = 30 mm 
(A = 0.39 and a = 2 ± 0.3). This indicates that the rear- 
rangement depends mostly on the normalized displace- 
ment z/D. The variation of M app with M also verifies 
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FIG. 6: Variation of the normalized apparent mass 
Afapp/Afapp with the normalized displacement z/D, 
at velocities 0.4mms _1 < v <34mms _1 (M = 300g, 
M^ = M 00 =98g). 

Eq. ©: the corresponding shielding length is indepen- 
dent of velocity with A = (48 ± 4) mm. As expected, 
A increases with D; experiments for a broader range of 
values of d and D will be necessary to establish the de- 
pendence of A on these parameters and on the ratio D/d. 

Preliminary experiments were performed at a lower 
relative humidity (H <40%). On the one hand, the vari- 
ations of M app with the total mass M remain in agree- 
ment with Janssen's model and both the mean M app val- 
ues and their relative dispersion are still much smaller 
than for the initial static packing. These mean values 
are also of the same order as before. On the other hand, 
M app remains constant in region (2) instead of increasing 
towards a limit value as previously. At the same time, 
grains remain at rest on the outside of the packing while 
the tube moves. These results imply that friction forces 
on the grains are smaller for H ~ 40%; they also indicate 
that even invisible rearrangements of the packing during 
the wall motion may induce large changes of the force 
distribution. Note that all experimental data presented 
here were obtained after an equilibrium of humidity be- 
tween air and the grains had been reached. 

These results indicate that the grain packing reaches 
a dynamical equilibrium independent of the initial state 
and of the relative velocity with respect to the walls: 
frequent transitions between different force distributions 
in the packing are induced by the tube motion leading 
to a well defined average value. These distributions are 
sensitive to small motions of the grains: these do not 
induce measurable variations of the mean particle frac- 
tion and are not even observable visually in some cases. 
The dynamical equilibrium is reached exponentially with 
a characteristic length proportional to the tube diame- 
ter: this may be due to the fact that the reorganization 
of the grains starts at the walls and propagates there- 
after towards the center of the tube. Stopping the flow 



increases M app , particularly at large velocities, but it re- 
mains smaller than the initial value after filling. The 
variability of M app from one experiment to another is 
much greater than during the tube motion since the sys- 
tem is "frozen" in one state and the averaging effect of 
the variations of the force distributions has disappeared. 

To conclude, the present experiments demonstrate 
that Janssen's model remains valid up to velocities of 
several centimeters per second for granular packing in- 
side a vertical tube in motion with respect to the grains. 
The apparent mass M app measured at the bottom of the 
packing during the tube motion is much smaller than for 
the initial static packing. The limiting value of the appar- 
ent mass at the end of the moving phase (M^p 1 ) is also 
independent of the tube velocity v and much less variable 
from a sample to another than for static packings (for a 
given total mass of the grains). Globally, the result of this 
study suggests that force distributions in dense granular 
flows may be described by straightforward models over a 
much broader range of velocities than usually expected. 
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